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ABSTRACT 

The overall purpose of this work is to develop a sub- 
sonic, turbulent, diffuse electric discharge operating at 
near atmospheric pressure for use as the pumping mechanism 
in an Electric Discharge Convection Laser (EDCL) . The 
work was done in two phases. The first consisted of a 
data collection phase during which the variation of dis- 
charge parameters was studied to determine their effects 
on the current-voltage characteristics of the multi-pin 
electrode arrangements used. The parameters studied 
included gap length, flow speed, turbulence intensity, 
electrode orientation, rate of voltage increase, and gas 
density through the addition of a diffuser. The second 
phase consisted of developing an interactive model which 
correctly reflects the dependences observed in phase one. 

A space-charge model resulted which correctly explains 
most major observed trends in gross discharge properties. 
Discharge stability is also discussed in terms of the model. 
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I . INTRODUCTION 



A. GENERAL 

This work is the continuation of an ongoing project at 
the Naval Postgraduate School (NPS) where the feasibility 
of constructing a high-energy, high pressure, continuous- 
wave Electric Discharge Convection Laser (EDCL) is being 
investigated. 

Laser systems of this t 5 rpe currently under investigation 
are limited in power output by the amount of electrical 
power which can be coupled into the laser (Ref. 1) . Since 
laser gain is directly proportional to the population 
inversion created by electrical pumping, the importance of 
optimum conditions in the discharge region can readily be 
seen. 

Previous work on this project at NPS vividly demonstrated 
the importance of flow conditioning on discharge performance. 
Power increases in excess of 250 times the normal no -flow 
conditions were demonstrated by flowing air rapidly through 
the discharge region with proper conditioning by turbulence 
generating screens (Ref. 2) . 

A brief general look at the physical situation pre- 
vailing in the discharge section will be presented here 
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to review the state of the art and to aid in interpretation 
of observed electrical characteristics. The basic geometric 
design of the electrodes is a common one often used in 

I* 

electric discharge work, the positive point-to-plane. This 
design results in an electric field distribution in the 
absence of space charge which is proportional to 1/x, where 
X is the distance from the point. This inverse field 
dependence is common to many non-uniform field electrode 
geometries including the hyperboloidal point-to-plane, 
concentric cylinders, wire-to-plane, sphere-to-plane, and 
many others, so the analysis presented here can be applied 
to a number of technically interesting cases. 

The field inverse dependence on distance from an elec- 
trode results in the division of the electrode gap into 
distinct regions of very high field where x is small, and 
very low field when x is relatively large. The smaller 
the radius of the electrode involved, the more pronounced 
the division. Figure (1) illustrates this effect for two 
geometries. 

As a result of the above, as voltage is increased across 
the gap, ionization will begin in the high field region, 
confining the discharge activity to the immediate vicinity 
of the pin where free electrons are able to gain enough 
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energy to ionize by collision. Just prior to complete 
breakdown of the gap, there exists a very small region of 
intense ionization which appears as a bright glow at the 
pin tip (in a point-to-plane discharge), and a space charge 
region of larger proportions, caused by less intense ioniza- 
tion and electrostatic repulsion of the positive ions formed 
in the ionization region by the positive pin. Beyond this 
space charge region lies a quasi-neutral region which con- 
sists of most of the low field region of Figure .1. The 
cathode, while lying well outside the high field region, 
does not play an important role in the dynamic processes 
taking place near the anode, and so will not be discussed 
in this work, except in its role at final breakdown. 

B. THE MODEL 

The analysis which follows is based on the ability 
of the flow field to alter the shape and extent of the 
space charge region by convection and turbulence. The 
problem is complex and involves many variables including: 

(1) tip radius (R) 

(2) gap length (d) 

(3) mean gas density ( p) and composition 

(4) mean flow velocity (U) 
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(5) electric/flow field coupling mechanisms (E,U,k) 

(6) time dependent fluctuations of 3 and 4 above 
due to the introduction of turbulence, U(t) , 

p(t) 

(7) ionization and charged particle djTiamics (a , 
a , X , e, n) 

The result is that the easily measurable gross discharge 
electrical properties such as voltage and current are 
functions of these variables, i.e«; 



I=f^(R, d,p , U, E, k, t) 
V=f2(R, d,P , U, E, k, t) 



( 1 ) 



There are many characteristic times associated with the 
variables cited, such as gas transit times, characteristic 
times associated with turbulent fluctuations of density 
and velocity, characteristic ionization and recombination 
times, streamer propagation and growth times, and many 
others (Table I) . The times cited above, however, will be 
of primary interest in the discussion which follows. 
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TABLE I 



Process 

Ionization 

Excitation 



Characteristic Time (Sec) Reference 
10 ’ 
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Energy Transfer from 
Electrons to Atoms -- 



Maxwellian Distrib. 
of Atoms 



Recombination 

Current Pulsing 

Streamer Propagation 
Turbulence Freq 
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Rate of Voltage Increase 
( V/y) 



-12 


10 


-12 


10 




- 10 


1 

00 

1 

o 

1 

0 

1 
1 


10 


■^-10"^ 


13 


■^-10"^ 


4 


-8 


4 


-1-10-5 


3 



0-2.5(10"^) 



(present 

work) 



(NOTE: The times from Ref. 10 are for an arc discharge.) 



Since the model presented in the following sections 
depends on the ability of the flow field to couple with the 
charged particles produced in the discharge region, it is 
important to understand how this coupling may take place. 
The three t 5 T>es of charged particles of primary interest 
are, of course, positive and negative ions, and free elec- 
trons, The bulk of the gas remains un-ionized over the 
major portion of the discharge. 
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Charged particles in the gas xmder the influence of 
an electric field acquire a drift velocity (V^=kE) which 
depends on the gas composition, pressure, and electric 
field strength. If the gas is flowing, then the particle 
velocity (in vector form) relative to the stationary 
electrodes is: 

"^total= kE + U (2) 

The sign convention used here will be that E, U and x are 
positive from anode to cathode. The k in Eq. (2) is the 
charged particle mobility, which can be expressed in terms 
of the parameters listed as k 4 e x/mv, where m is the 
particle mass, X (the particle mean free path) is inversely 
proportional to gas pressure, and v is the random thermal 
velocity of the particle. Equation (2) can be rewritten 
as = kE(l + U/kE) , where the quantity U/kE will be 

referred to as the "slip parameter, S.” Electrons, due to 
their small mass, have a very high mobility, and generally, 
the slip parameter for electrons is very small, and gas 
convection is not significant; but ions are much more 
massive, so k is low and ion drift velocities can be of 
the same order of magnitude as U, even in the high field 
regions of the discharge, resulting in significant ion 
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relocation by convection within the characteristic times 
mentioned above. Since high densities of charged particles 
produce electrostatic fields of their own, which can easily 
reach the same order of magnitude as the applied field, the 
electrons are affected by convection indirectly through 
the alteration of the electric field structure by the move- 
ment of high density ionic space charges. 

Visual observations of a discharge vividly demonstrate 
the ability of highly turbulent flow to affect the space 
charge. With high velocity and no turbulence, the glow 
extending from the tips of the pins toward the cathode is 
very bright and steady, while, when high turbulence is 
introduced, the bright portion of the glow can be seen to 
radically "whip" around, indicating the high degree of 
influence that the fluctuating velocity components have 
on the charges present in the flow. 
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II. EXPERIMENTAL APPARATUS 



The primary test apparatus used in this work consists 
of a plexiglass wind tunnel with a test section in which 
various electrodes and turbulence generating grids were 
mounted. Electrode spacing is adjustable and the electrode 
orientation with respect to the flow field could be changed, 
permitting much flexibility in the examination of flow and 
electric field interactions for a variety of conditions and 
with several levels of turbulence. Later in the testing, 
a new wind tunnel was added which was equipped with a 
diffuser to permit operation at higher flow velocities and 
reduced density in the discharge section (Ref. 3) . 

The electrodes are of a design common to many electric- 
discharge lasers (Fig. 2) . The anode consists of a pin 
rack with three rows of pins, eleven to thirteen pins per 
row. The pins were made of steel, with horizontal and 
vertical separation of one centimeter. The cathodes used 
varied, because there was a need to search for the most 
effective and structurally soiond arrangement. All the 
cathodes used, however, were geometrically similar, and 
simulated a plane with minimiam solid area to reduce flow 
blockage as much as possible. Structural integrity became 
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an important consideration as testing progressed, since 
flow velocities exceeded 180 meters per second and turbulence 
levels introduced were relatively high (Ref, 3). 

All tests were performed in an air flow which was 
supplied by a Carrier three-stage Centrifugal Compressor 
with a maximum inlet capacity of 4000 ft /min, and a maximum 
pressure ratio of two. Flow velocity was regulated through 
a series of valves, and flow total temperature was kept 
below 100°F by passing the air through an air cooler prior 
to reaching the test section. Flow velocity was measured 
with pitot or pitot- static tubes which were mounted in or 
downstream of the test section, depending on the tunnel 
being used. Turbulence was introduced by movmting one or 
more grids upstream of the discharge section. Some grids 
which were used in previous work were retested for verifi- 
cation and validation of systems operation (Ref. 3), and 
several new designs were tested (Fig. 3) . 

Turbulence spectra were measured with a Thermo Systems, 
Inc. hot-wire anemometer, which consisted of a 1051-2 monitor 
and power supply, 1054A-30 anemometer module, and locally 
manufactured hot-wire probes, A detailed account of turbu- 
lence measurement equipment and techniques can be found in 
Reference (3) . 
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Two power supplies were used to power the electric dis- 
charge. The first was a Sorensen Model 1030-20 High Voltage 
D.C. power supply, which is rated at 30 kV DoC. and 20 ma 
with a 2% ripple. It was soon discovered that the maximvnn 
gap length was severely limited by this supply, so a Sorensen 
BETA Model High Voltage Supply was obtained which was rated 
at 60 kV D.C. and 50 ma (Fig. 4) . With this power supply 
it was hoped to increase the maximimi gap length which could 
undergo breakdown. Due to \mcorrectible internal problems, 
the maximum voltage and current obtainable was 34 kV and 19 ma, 
so reliable tests were limited to less than 4 cm gap length. 
Both power supplies were equipped with trip controls which 
work on either voltage or current. The voltage was initially 
measured by a Sensitive Research electrostatic voltmeter 
with a maximum range of 40 kV and an internal impedance of 
5.10^^ ohms. The readings of the power supply voltmeters 
were compared and were found to be in excellent agreement, 
so they were used exclusively throughout the remaining tests. 
The current was measured with a variety of ammeters, ranging 
from microamp scales to a 200 ma scale. A meter panel was 
set up with three ammeters in parallel, one 0-1 ma range, 
one 0-50 ma, and one 0-200 ma. The meters were wired 
through a multi-position "Make-before-break” switch so that 
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as the range of one meter was exceeded, the next one could 
be switched in without circuit interruption. A Simpson 
volt-ammeter was also wired in series with the ammeter 
panel for use in the microampere range and as a check against 
the readings of the panel meters. This redundancy permitted 
immediate detection of erroneous current readings, and re- 
placement of defective meters. All wiring was done with 
high voltage leads which were fabricated prior to the initial 
test phase. Connections were made through highly polished 
brass spheres which were about 1% inches in diameter. The 
high voltage side from the power supply outlet to the anode 
was immersed in a high dielectric oil bath to prevent 
voltage leaks at the ball connection. Great care was taken 
to ensure that all equipment was well grounded at all times, 
through the laboratory grounding system. 

One series of tests performed required increasing the 
applied voltage across the gap at various rates. A constant- 
speed 24 volt D.C. motor, powered by an adjustable output 
D.C, Power Supply was connected to the Sorensen 30 kV 
power supply voltage control shaft through a Matron Variable- 
Speed transmission (Fig, 5) . This arrangement permitted the 
applied voltage to be increased across the gap at rates from 
0-20,000 volts/second. 
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III. EXPERIMENTAL OBJECTIVES AND MODELING APPROACH 



A. GENERAL 

A brief preview of the experimental work and an overall 
perspective of the model development should serve to help 
keep sight of the "big picture" while sorting out the de- 
tails of the individual pieces. 

The experimental work has been divided into two main 
categories. First, tests were run varying individual dis- 
charge parameters, including gap length, flow velocity, gas 
density, turbulence intensity, cathode configuration, and 
electrode orientation. The purpose of these tests was to 
establish a data base for use in analyzing the individual 
contributions of various parameters to discharge power, 
and to attempt to formulate an interactive physical model 
based on the observed relations. The second phase of test- 
ing consisted of observing the effects of increasing applied 

voltage across the discharge gap at various rates to deter- 
«. 

mine what effect this would have on the breakdown properties 
of the gap. The rate of voltage increase was constant for 
each run, but it varied from run to run. These tests also 
added to the data base which was being established, but the 
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specific purpose was to determine if any type of resonance 
effect might be detected between a particular rate of in- 
crease, and the turbulence spectra which were being intro- 
duced since the characteristic time scales of turbulence 
and the rates of increase were of the same order of magnitude. 

B. MODELING APPROACH 

A space-charge model development, which begins in 
Section IV, is based on the dynamic processes which occur 
when the gap is very near its breakdown potential. Under 
these conditions, a positive point-to-plane gap in an 
electronegative gas such as air, experiences a pulsating 
current due to the formation of pre-breakdown streamers 
whose magnitudes vary in a statistical way. The frequency 
of these pulses has been measured on the order of 10^ to 
10^ Hz (Ref. 4) . The model proposed here looks at the 
discharge from the point of view of an individual event, 
such as the propagation of a pre-breakdown streamer, and 
the situation that prevails in the discharge at the time 
and location of its occurence. For example, in a steady, 
non-flowing point-plane discharge near breakdown potential, 
the streamers which cyclically form, all "see” generally 
the same "picture" (field structure, space charge distribu- 
tion, etc.) as the preceding ones did. Much activity takes 
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place between these events, but at each period the situa- 
tion at a given location is very similar. Now, if some- 
thing is done to alter the picture that the streamers see, 
such as a reduction of density, or flowing gas through the 
gap, or the introduction of turbulence, the "picture" pre- 
sented to each streamer may be very different from the 
steady, no-flow case. The change in the local situation 
can be the result of two separate classes of influencing 
factors, one steady-state, such as a permanent change in 
gas density; and the other, a time- dependent change, such 
as with the introduction of convection or turbulence, where 
the total change seen by the next streamer depends on how 
much time the effect has to work on the discharge between 
streamers. The change in this picture over the original 
case, then, is what the model assumes determines the changes 
observed in gross discharge properties. Of course, both 
classes of influence can be operative at once, and then 
their individual effects must be sorted out. 

To simmarize, the model was developed by isolating a 
particular influence factor, determining its effect on the 
"picture" seen by successive streamers, then attempting to 
combine the results of all such influences in an interactive 
way to explain observed changes in gross discharge 
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characteristics. The ultimate goal of this study is to 
develop a method of approach for optimizing the discharge 
for use as a pumping mechanism for the EDCL. 

C. PREVIEW OF EXPERIMENTAL/MODEL CORRELATIONS 

The first part of Section IV introduces the basic space- 
charge model and illustrates how a reduction in density by 
addition of a diffuser (Ref. 3) alters the picture presented 
to succeeding streamers. Successful development of a simple 
relation for the density dependence of gap breakdown voltage 
in a non-uniform field evolves from fitting existing data 
to an equation form resulting from the space-charge model. 
This lends credibility to the basic concept. 

The next part shows an analysis of the influence of 
convection on the discharge with various orientations of 
electric and flow fields. With the aid of the space-charge 
model again, the change in "picture" seen by successive 
streamers is determined, and various maxima and minima for 
breakdown voltages are calculated for variations of velocity 
and orientation. Again the results obtained agree with 
experimentally observed values, lending further credibility 
to the proposed model concept. Various instabilities and 
observed visual characteristics are also discussed in terms 
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of the model, with general qualitative agreement between 
model predictions and observations. 

In the second half of the development, velocity, density, 
and turbulence effects on discharge current and the inter- 
relation between the current and voltage are discussed. A 
current equation evolves from the space-charge model which 
properly accoxmts for most major observed trends in dis- 
charge current due to convection, turbulence, and density 
reduction in various flow orientations. Some stability 
conditions are mentioned which are particularly relevant 
to the discussion, but no attempt is made to give a detailed 
account of stability effects. 

The four basic discharge conditions used for comparison 
of the space-charge model predictions are the following; 

(1) No diffuser, no turbulence 

(2) No diffuser, with turbulence 

(3) Diffuser on, no turbulence 

(4) Diffuser on, with turbulence 

Test results appear to bear out the predicted trends for 
each of the above cases. Detailed design data on the dis- 
charge apparatus, including diffuser and turbulence 
generators are given in Reference (3) . 
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IV. SPACE CHARGE MODEL 



Ntnnerous studies have been done on various aspects of 
flowing gaseous discharges in order to determine the effects 
of convection and/or turbulence on the strength and stability 
of such discharges (Refs. 2, 3, 5, 6) . 

As pointed out in the Introduction, xjnders tending the 
basic processes in a non-uniform, flowing discharge is 
essential to the improvement of the ultimate potential of 
the Electric Discharge Convection Laser. Present discharge 
stability problems restrict practical operation of this 
type of laser to the realm of 100-200 torr operating pres- 
sure (Ref. 5) . This pressure limitation, or more accurately, 
density limitation, in turn limits maximtmi power input to 
the gas, and, consequently, restricts laser power output. 

In the following discussion, an attempt is made to shed 
some light on the complex problem of discharge stability 
and power enhancement in a non-uniform, subsonic, near- 
atmospheric flowing discharge. Using some well-known 
principles of gaseous breakdown in uniform fields and a 
few basic gas dynamic ideas, it is possible to, at least 
qualitatively, explain many diverse observed phenomena of 
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the particular discharge sections used in this research. 

The basic concepts, however, may be extended to any geometry 
and any flow condition, provided of course that these differ 
ences are taken into account in interpreting observed dis- 
charge characteristics. 

A. BREAKDOWN VOLTAGE (V, ) 

D 

1. Uniform Field (parallel plate) 

Cobine (Ref. 7) and von Engle (Ref. 8) both analyze 
the effects of space charge on the current and voltage dis- 
tributions in the gap between parallel plate electrodes. 

The case of the positive ion emitting anode from Reference 
(8) is suitable for the purposes of the present discussion. 
Consider first a uniform electric field between the plates 
with no space charge. The electric field will be constant 
from anode to cathode at: 

E=V/d 

where V=applied voltage 

and d=electrode spacing. 

Now, consider the effect of a space charge modifica- 
tion of the field structure, with the anode emitting posi- 
tive ions. Assuming that the bulk of the discharge is 
electrically neutral as is normally the case, there will 
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be a small region near the anode in which the ion density 
exceeds the neutral charge density. A net positive charge 
will then characterize this area. 

The effect of the positive space charge is to lower 
the electric field in this region, resulting in the type 
of field and voltage distributions shown in Figure 6. The 
modified E and V curves assxme constant charge density out 
to X=Xg^, which is not physically possible, but it illus- 
trates the point. It should be noted from the figure that 
an intense space charge can actually cause a positive 
initial voltage gradient. At point P in the figure, the 
voltage gradient is zero, and point P is referred to as a 
"virtual" anode. 

Referring to the distance at which the virtual anode 
lies from the origin as 6 , it can be seen that the value 
of the electric field in the neutral region must increase 
to account for the same voltage drop between the electrodes. 
One way of interpreting this is to define an equivalent gap 
length 

d'=d-6 (3) 

Noting also that the effective voltage drop from the 
virtual anode to the cathode is increased by the amount AV, 
it is possible to write an expression for the effective 
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electric field strength in the neutral region as: 

E« = V+ AV 
d' 

The important point here is that the apparent physical gap 
length is reduced by space charge modification of the 
electric field structure (Ref, 9) , This "virtual" electrode 
effect applies to any electrode near which charges of the 
same polarity accumulate (Ref. 7) . 

In a normal parallel plate gap, the situation is 
very similar to the one just described, since positive 
charge tends to build up around the anode (Ref, 4) . If 
this charge reaches a high enough concentration just prior 
to breakdown, the result will be a "virtual" anode, and the 
effective gap length d‘ will reflect this reduction. 

Neglecting AV, then, the breakdown field strength 
would be 

% ^ V^/d’ = V^/(d-6 ) 
which can be written as 

Eb = (Vb/d)/(l- Vd) 

and for <5/d<< 1, this can be written as: 

= (V^/d)(l + 6/d) (4) 

This equation provides a basis for a test of the space-charge 
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model concept which will be used to describe processes in 
the non-vmiform field, so it is worthwhile to see if it 
gives results approximately in agreement with observations. 

Since it was assumed that 6/d was small, Ej^ can 
be evaluated at as large a value of d as possible to mini- 
mize the contribution of the <S/d term. From the data of 
Reference (4) , at d=20 cm, V^=510 kV, so 

= \/d)^ = 510/20 = 25.5 kV/cm 

Now, to estimate <S, the data point at d=10 cm can be used, 
where the measured V^=265 kV. The hypothesis is that for 
any gap, 

Eb = V^/d)^ (1 + 6/d = 265/10 

so, for the 10 cm gap, 

(1 + 6/10) = 26.5/25.5 

and 

6 = .39 cm 

Now, the final result is: 

E^ = + 6/d) = 25.5 (1 + .39/d) (5) 

In reality, the anode extension can never be larger than 
the physical gap, so the upper limit of Ej^ based on this 
crude approximation, is 2 E^q^ = 51kV. 
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It is also noted that in the absence of space charge 
effects (i.e. 6 =0), E^=const which is the required result. 

It can be seen in Figure (7) that at 2 cm and above, 
the computed values of the breakdown field are quite close 
to the measured values; closer, in fact, than the values 
computed from the other models listed in Reference (4) . 

This agreement demonstrates the possible utility of such a 
space-charge model, and this concept will be used in a 
detailed analysis of gasdynamic effects on a non-uniform 
field. 

2 . Non-Uniform Fields 

The basic concept described above for the uniform 
field can also be applied to the non-uniform field distri- 
bution used in the present experimental investigation. 

First, the process will be described as in a non-flowing 
discharge section, and later, the effects of gas convection 
and turbulence will be discussed. Experimental data will 
then be scrutinized in light of the arguments presented. 

a. No -Flow 

First, consider a stationary discharge in a 
point-plane gap. For sharp points, the region near the 
point is characterized by very high electric fields just 
prior to spark breakdown of the gap. For a simple model. 
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consider the field distribution for a hyperboloidal point- 
to -plane (Ref. 10) : 



E(x) 





2V 

R log(4 d/R) 



( 6 ) 



Ep = field at surface of point 
R = point radius 
V = applied voltage 
d = gap length 

For typical values of R and d of 0.05 and 3 cm 
respectively, breakdown voltage has been found experimentally 
to be about 27,400 volts in atmospheric air. This results 
in an Ep of: 

E = = 200,000 v/cm 

P 0.05 log(4-3/.05) 

This corresponds to a reduced field (Ep/P) of 263 v/cm-torr. 

Since this value of electric field is well above 
the atmospheric air ionization field strength of 15,000 
v/cm (20 v/cm-torr) (Ref. 11) , it is evident that intense 
ionization of the air will be taking place in the vicinity 
of the point by high energy electrons accelerating toward 
the positive point. 
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As an electron accelerates toward the point, 
it creates an electron avalanche. Since the electron 
mobility is much higher than the ion mobility, the elec- 
trons speed into the positive point and get absorbed 
leaving behind the ion space charge in front of the pin 
tip. Detailed descriptions of the development of streamers 
and the electron avalanche can be found in References (7, 

8, 10, 12). It is sufficient for the present to realize 
that dense positive space charge can be formed in the 
immediate vicinity of the anode which can significantly 
alter the electric field structure there. (In an electro- 
negative gas such as air, negative ions would also be 
formed, and their number, to a first approximation, can be 
regarded as a fraction of the number of positive ions.) 

The net charge density in the vicinity of the anode, there- 
fore, remains positive, and negative ion effects will be 
disregarded for the time being (Ref. 9). The positive ions, 
under the influence of electrostatic repulsion and diffusion, 
begin to move slowly away from the positive point, creating 
the space charge region. As voltage is raised, the dis- 
charge passes through a stable, pulseless "glow" mode, and 
then, as applied voltage approaches the gap breakdown value, 
the formation of pre-breakdown streamers causes the frequency 
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of current pulses to increase to approximately 10 to 
10^ Hz. The maximum length of streamers measured under 
pulsed conditions has been observed to be a linear ftmction 
of the applied voltage when no space charge exists to dis- 
tort the electric field (Ref. 4), It seems reasonable, 
then, that the pre-breakdown streamers formed in a static, 
steady-state discharge would behave in a similar manner. 

The high concentrations of positive charge in the vicinity 
of the anode cause an apparent anode extension into the gap, 
reducing the effective gap length much as in the case of the 
’’virtual" anode described in the previous section. These 
pre-breakdown streamers have only to travel the reduced, 
effective gap, then, to effect breakdown of the gap by the 
processes described in Reference (4) . If this assumption 
is correct, then influencing the extent of the space-charge 
region by flow or turbulence should show up as changes in 
the measured values of breakdown voltage. To test this 
assxnnption, it is necessary to estimate what the original 
size of the space charge region is \ander no-flow conditions. 
The following sections address themselves to this problem, 
and then to the problem of determining the changes in 
breakdown voltage which might be expected due to alteration 
of the space-charge region by gasd 3 mamic processes. 
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(1) Spatial Development of the Virtual Anode , 



It is interesting, and will be useful later, to examine 
the effects of gas density variation on the anode exten- 
sion and breakdown voltage under no-flow conditions. The 
following relationships are used: 



Mean Electron Energy 
for Ionization: 


eE = eV^ 


(7) 


Electric Field 
Distribution: 


E(x) = 


(8) 


Breakdown Voltage 
Relation: 


V^ ± ad'+b 


(9) 


Virtual Anode Effect: 


1 

II 


(10) 


Mean Free Path: 


^ = /.\S~ \ ^ c 


(11) 


3 P(torr) 



where 

-19 

e = electronic charge = 1.6x10 coul 

= ionization potential = 15.5 volts for air 
a,b,c = constants 

X = electron and gas molecule mean free paths 
® »8 

P = pressure 

and all other symbols have been previously defined. 

Solving for the ionization field strength 

from Eq. (7) , 

Ei = V./ X^ 

> 
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and using Eq. (11), this becomes 

= V^P/c 

Now, using this value in Eq. (8) , it is possible to solve 
for a characteristic "ionization distance" (the maximtmi 
distance from the anode at which a free electron would, on 
the average, in one mean free path, gain sufficient energy 
to ionize by collision). Rearranging Eq. (8) with E(x) = E^, 
and recalling E^ from Eq. (6): 



(R+x^) 



R p = R Ep(v) 4 jv 
Ei P (V^P/c) P 



It is assumed that the anode extension 6 is proportional 
to x^, and that R<<6 , so that 



(S i X. = JL = JL. (Const temp) 
i P p (Ideal gas) 



( 12 ) 



Now, in direct analogy with the uniform field example of 
the previous section, the effective gap length is expressed 
by Eq. (10) as 

d' = d-6 d- — c' 

P 



and it is seen that the effective gap length is reduced as 

the gas density is reduced. If this value of d' is inserted 

into Eq. (9), after some rearrangement it can be shown that 

V = (ad+b) p /i2) 

^ ( (H- ac') 
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The form of this equation will be shown to describe the 
actual breakdown voltage density dependence fairly well. 
For example, if it is rewritten as 



where 






P / Pn) 

( p/ Po+1) 



(14) 



Vbo = 



at one atmosphere 



the values of vs. P / Pq agree with the observation. 

The constant term in Eq. (14) is dependent on the specific 
electrode geometry and type of gas. Table II is a collec- 
tion of data for a positive point -to -plane for gap lengths 
and densities in the present range of interest. 

In making comparisons of measured data at 
one atmosphere pressure, it was noted that Bandel's data 
from Table II for the 8 cm gap is inconsistent with the 
data of References (6) and (7) for a point-plane gap of 
8 cm, and with trends noted in the present experimental work. 
Also, due to the fact that Bandel's breakdown voltage at one 
atmosphere was nearly the same for 4 and 8 cm gaps, there 
is doubt that the value of given is appropriate for the 
present analysis. Therefore, the value of used in the 
calculation of the voltage ratio for the 8 cm gap for com- 
parison with Eq, (14) was 60 kV, which is in close agreement 
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with the other references cited. Bandel's second point 
for the 8 cm gap was also discarded for the reasons already 
mentioned. Figure (8) shows the results of observed values 
of the voltage ratio vs. the predicted values of Eq. (14). 
The close agreement tends to support the "virtual" anode/ 
space-charge effects model outlined above. 



TABLE II 



Density 

Ratio 


Vb 


(kV) 


Ib ( 


a) 








Meas 


/ Calc 


Meas / 


Calc 


Gap 

Length 


Source 


.05 


2.0 


2.6 


0 


0 


4 


Bandel 


.09 


4.0 


5.8 


65 


60 


4 


(Ref 18) 

It 


.15 


7.2 


8.7 


105 


110 


4 


It 


.24 


14.2 


14.7 


225 


180 


4 


It 


.39 


22 


21 


240 


230 


4 


ft 


.63 


33 


33 


240 


230 


4 


It 


.98 


38 


38 


130 


130 


4 


ft 


.05 


4.0 


5.4 


25 


20 


8 


11 


.09 


8.0 


9.0 


150 


90 


8 


It 


.15 


15 


14.3 


180 


166 


8 


It 


.24 


24 


21.8 


215 


204 


8 


ft 


.39 


32 


31 


170 


185 


8 


II 
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TABLE II (Continued) 



Density 

Ratio 


Vb 


(kV) 


Ib 


( a) 




• 




Meas 


/ Calc 


Meas 


/ Calc 


Gap 

Length 


Source 


.63 


(33) 


40 


105 


105 


8 


Bandel 


.98 


(38) 


41 


40 


36 


8 


(Ref 18) 

It 


1.0 


4.8 


2.2 


8.5 


1.7 


.3 


Werner 


1.0 


5.0 


5.8 


5.5 


4.8 


.6 


(Ref 15) 

II 


1.0 


9.2 


11 


14 


10.6 


1.0 


Townsend 


1.0 


10 


18 


(10) 


21.5 


1.5 


(Ref 17) 

tf 


1.0 


25 


25 




37.5 


2.0 


Meek & Graggs 


1.0 


60 


-- 


— 


— 


8.0 


(Ref 10) 

It 


1.0 


65 








8.0 


von Engle 










(Ref 7) 



(See Appendix B for calculation procedures and formulae.) 

(NOTE: Data in parenthesis do not seem to correlate with 

trends observed in other references.) 



b. Gas Convection/Flow Orientation Effects Analysis 
It has been observed from our work that the 
orientation of the flow field relative to the electric field 
has a marked effect on the electrical properties of the 
discharge. Most evident is the change of breakdown voltage 
observed in the various orientations tested (Fig. 9). 
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Figures (10, 11) illustrate the nature of these changes. 

It can be seen from the Figures that convection increases 
Vb/Vbo 111 the reverse and cross-flow modes (Fig. 11) 
quickly to a nearly constant value on the order of 1.25. 

In the ordinary mode, however, the general trend is for an 
initial decrease in to a minimum in the neighborhood 

of 0.9-0.95 at aroxand 75 m/sec flow velocity. After that, 
the ratio begins to increase again, occasionally to a 
value greater than one. 

A space-charge-convection-effects model has 
been developed which successfully explains some of the 
major features of the breakdown voltage behavior. 

Two orientations were tested, one longitudinal 
mode with flow from anode to cathode (ordinary mode) (Fig. 
9a) , and one cross-flow mode with the discharge being 
transverse to the gas flow (Fig. 9c). The third mode 
analyzed by the present model is another longitudinal mode, 
as in case one, but with the flow direction reversed so 
that the gas flow is from cathode to anode (Fig. 9b) . 

(1) Space-Charge-Convection-Effects Model . 
Although the discharge sections used consisted of from 33-*^ 
39 positive pins spaced one centimeter apart for the anode, 
and a grounded plane- t 5 rpe cathode, the electric field 
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structure for any given pin is nearly identical with that 
of a single point-plane configuration, since the pins are 
far enough apart that the space charge regions do not over- 
lap (Fig. *2), Therefore, the model is discussed in terms 
of the d 3 mamic processes taking place at one pin, with the 
understanding that the processes apply equally to the test 
section as a whole, provided that all pins contribute 
uniformly to the current. 

It must be noted at this point, however, 
that due to the statistical nature of breakdown streamer 
size and intensity, breakdown becomes a probability related 
process. For example, if the voltage across a single point 
gap is set so as to give a probability of breakdown 
in a certain time interval (At) of Pf]^~ol, then its prob- 
ability of not breaking down is 

(15) 

Now, if there are 10 pins acting, inde- 
pendently at the same applied voltage, each with 
then the probability of the 10 pin gap not breaking down 
would be given by the following expression: 

Psio= l-(lPfi)’-° - = .65 (16) 

This demonstrates the increased susceptibility of a multi-pin 
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configuration (with statistically independent pins) to 
break down. Therefore, to obtain the same breakdown prob- 
ability, a multi-pin gap must be at a lower voltage, and 
observed values of breakdown voltage in multiple pin con- 
figurations may reflect this reduction. 

The expression used for the electric field 
in the space-charge convection model is the same as used 
in previous sections i.e., that for the hyperboloidal point- 
to-plane. In this model the positive space charge is 
assumed to be at a high enough density to cause significant 
distortion of the undisturbed field out to a distance where 
the undistorted positive point field strength is - 20,000 
volts/cm. This point in the field was chosen because it 
is well inside the region of ionization by electron collision 
of = 15,000 volts /cm for air at atmospheric pressure. The 
corresponding reduced field strength (E/P) = 20 v/cm-torr. 
Referring again to Figure (1) , an electron arriving at point 
a, E = 15,000 v/cm, would, on the average, have enough 
energy to ionize air molecules by collision and electron 
avalanching would start toward the positive tip. The 
result, as described in the previous section, is positive 
space charge acctmulation in front of the pin. The posi- 
tion where the pin field strength approaches 20,000 v/cm. 
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point b, is conjectured to be in the critical region where 
the space charge field becomes significant compared with 
the undistorted field. When this happens, the rest of the 
gap acts as if the anode lies inside the high field region 
caused by the dense positive space charge, and the gap is 
effectively shortened. Figure (12) is an illustration of 
this idea. 

The distance defined above is used in the 
rest of the analysis as the magnitude of the anode extension, 
previously defined as 6 in the no-flow analysis of the 
density dependence of the breakdown voltage. The effective 
gap length, as proposed by the model, would be the distance 
between the cathode and the point at which the static elec- 
tric field first reaches 20,000 v/cm. Expressed mathematically, 
this would be 

d' = d - 

where the subscript indicates the field value at which 
is computed. This value of 6 will represent the no-flow 
anode extension and will be referred to as 6^= 6^0^ . 

(2) Experimental Inputs and Basic Model . 
Experimental data over a range of gap lengths \mder no-flow 
conditions at atmospheric pressure have yielded an approxi- 
mately linear relation between gap length and breakdown 
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voltage, closely approximating the point-to-plane curves 
of References (8, 10) (Fig, 13). 



The apparent reduction in gap length, . , 
can again be calculated from Eq. (8), using E(x) = 20,000 
v/cm, as 




^20K 



- R E p . -R 

20,000 



(17) 



but note that 



Ep 



2V 

R log(4d/R) 



Using the experimentally derived relation for breakdown 
voltage (Fig, 13) : 

= 7800d+4000 (volts) (18) 

With an estimated R=0.05 cm, then 

.05( 2(7800d + 4000) ) q9) 

6 (d) = »03 1p K(Ad/l^ ..05 

20,000 

TABLE III 



d(cm) 5 (cm) (S /d 



2 


.34 


• 


4 


.56 


.14 



The quantity 'S/d in Table 3 is a measure 
of the relative reduction in gap length for two gaps in 
the range of interest. It can be seen that the relative 
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reduction due to the presence of space charge over the 
2 cm difference is only about 37o. What this means is that 
the linearity suggested in the data may not be between 
and d, the actual physical gap length, but may be a linearity 
between and d' , the apparent gap as speculated earlier. 

Due to the small relative reduction of gap length, it is 
not possible to determine from the data alone which is the 
actual case. 

In the ensuing discussion, it is assumed 

that V, is a linear function of the effective, reduced gap; 
b 

= d' = d - 6 (9, 10) 

Adoption of this hypothesis leads to some interesting 
conclusions . 

(3) Calculation of Maximum Voltage Ratio due 
to Convection . In flow orientations two and three, (Fig. 

9 b&c) , it can be deduced from the above theory, that the 
maximvim value of breakdown voltage occurs when convection 
effects remove all the space charge from the gap, thereby 
increasing the effective gap spacing from its normal no- 
flow value of d^ = (d - 6^) to d, its upper limit. The 
max im u m ratio attainable then at any velocity should be 

= - t ~ r 

o 

where the subscript 'o’ denotes the no-flow value = )• 

O JLO 
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The Table below summarizes the results of 
this calculation for several gaps in the range of the ex- 
perimental data. 



TABLE IV 

REVERSE AND CROSS -FLOW MAXIMUM VOLTAGE RATIOS 

Gap Length Anode Ext Effective Gap Voltage Ratio (max) 

(Eq. 17) 

d(cm) 6 n ( cm) ^o Theory Experiment 



1 


CM 

CM 

• 


.78 


1.28 


-- 


1.4 


.27 


1.13 


1.24 


1.27 


1.45 


.27 


1.18 


1.23 


1.26 


2 


.34 


1.66 


1.20 




3 


.45 


2.55 


1.18 


1.20 


3.2 


.47 


2.73 


1.17 


1.15 


4 


.56 


3.44 


1.16 


1.13 



(See also Figure 14) 



It is evident from* the Table that theory 
and experiment agree quite well, despite many crude approxi- 
mations which went into the calculations. The close numeri- 
cal agreement tends to verify the initial assimiption that 
the point where E(x)=20,000 v/cm may be considered to be 
the approximate axial range of the anode extension under 
no-flow conditions. Perhaps even more significant than 
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the numerical agreement is the fact that the theory pre- 
dicts a decrease in the maximxjm voltage ratio with increas- 
ing gap length. This trend has been experimentally verified. 

(4) Minimxjm Voltage Ratio . The ordinary flow 
mode, Figure (9a) , is more complicated as can be seen from 
the data summary plot of Figure (10) . 

It might be expected from a space charge 
convection point of view, that in this orientation the 
space charge would be carried toward the cathode in some 
characteristic time interval which would tend to increase 
6 , If 6 increases, the effective gap length decreases, 
which in turn results in a decrease in V^, since - d'. 

This agrees with the initial trend of the data of Figure 
(10). The breakdown voltage vs. velocity curves show a 
minimum though, beyond which begins to increase again. 

At short gaps, in fact, has been observed to not decrease 
at all, but show a slight increase which remains relatively 
constant with increasing velocity. 

In order to shed some light on these 
apparent inconsistencies at short gaps, and to attempt to 
explain the complex nature of the breakdown in this orienta- 
tion, it is necessary to look into the relative magnitudes 
of the convective and field effects on the ions with 
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reference to their specific locations in the total field 
structure. A convenient parameter for such a comparison 
is the "slip parameter" 



g _ JU_ _ U_ _ gas velocity = convection effect (21) 
kE Vjj drift velocity electrostatic effect 

When the order of magnitude of the slip 
o 

parameter is = 1, convection effects become as important 
to ion motion and location as the E field, whereas if 
S = U/v^<< 1, convection effects would be negligible. 
Conversely, for S » 1, the electrostatic effects become 
negligible. Figures (15 & 16) illustrate the magnitude 
of this parameter plotted against U for the field point 
E=20,000 v/cm, and E=Ep, 

It is evident from Figure (15) that an 
ion located in an electric field of strength E=20,000 v/cm, 
would experience significant convective effects at very 
low velocities. Figure (16) is an illustration of the 
slip parameter values at the anode surface, S=U/v, , with- 

dp 

out considering the effects of space charge. However, it 
must be remembered that in actuality, a positive ion at 
the anode surface is located between the tip and a region 
of high positive space charge, which reduces the effective 
field strength at its location (See Ep' in Fig. 12). 
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By reducing the effective value of Ep by 

amount f (n**^ , then E '=E -f(n’*’) , and S =U/k'^E ' increases, 
* P P P P ’ 

which could render convection a significant effect even at 
the tips of the pins. An important aspect of the quantity 
U/k'^Ep* is that it is dependent on the gap length, since, 
as stated earlier, at breakdown: 



= 7800d + 4000 (volts) and 

E*, = - f(n"^) which yields 

Pb Rlog(4d/R) ^ 



_U U 

v^p’ k‘^2(7800d+4000) /Rlog(4d/R) 



Due to the nature of the log term, this yields approximately 
the following important dependence: 



S = 



i U 

Vdp' d-f(n-^) 






What this equation effectively says is that 
at short gaps, the convective effects will become more 
prominent and may tend to dominate the discharge, while at 
longer gaps, at least throv:gh the range tested, both factors 
may be important. It is noteworthy that the quantity f(n"*^ 
is proportional to the current, since the charged particle 
density is proportional to the current. This factor would 
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also tend to increase the slip parameter as the current 
increased. 

The alternative to attempting to analyse 
the exact magnitude of the effective electric field at the 
anode surface is to inspect the available data for evidence 
of significant flow effects which may alter the breakdown 
voltage characteristics due to creation of a sufficiently 
high slip parameter at the anode surface. 

A closer inspection of the data of Figure 
(10) indicates a reversal of the declining trend in vs. U 
for those gaps greater than about 2 cm in the range of 60- 
90 m/sec flow velocity. Therefore, it is assumed that 
this is the velocity region in which convection begins to 
seriously affect the ion population closest to the pins. 

The current vs. velocity curve also indicates that the 
current starts to rise sharply in this range, so the com- 
bined effects of high velocity, high current, and short 
gap (Eq. 22), appear to significantly alter the d 5 mamics of 
the situation. Since convection becomes important, it is 
reasoned that, in a certain characteristic time interval 
(The time between current pulses or breakdown streamers 
seems to be appropriate in this case, and can be estimated 
from the pulse frequency cited in the introduction as 
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'f: 




t'.’ 










4 - s 

-10 sec.)> the ions are displaced by such 



an amount away from the pin tip anodes, that the space 
charge no longer acts as an anode extension. Instead it 
becomes a detached region of floating charge, no longer 
able to act as a bridge for infeeding electrons at the 
head of succeeding avalanches, to run all the way into the 
anode itself. Figure (17) illustrates the various stages 
of progression of the events just described, A separate 
shortening effect of the anode extension, due to high 
current density and recombination, may also play a part 
in the observed voltage characteristic rise, and it will 
be discussed fully in the next section in connection with 
observations in the cross -flow mode. 



ordinary flow mode, the sequence of events as described is 
as follows: 



The physical gap is effectively reduced 
by the concentration of space charge. The effective gap 



To summarize up to this point for the 



(a) No Flow 



S = U/v^ = 0 



Fig. (17a) 



is d' = d-5 , and V, /V, =1.0. 

O O b DO 



(b) Intermediate Flow (IF) : 



® ' “/''d(20K)° ^ 



Fig. (17b) 
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Convective effects extend 5 further into , 



the gap and d^^, = d^ -UT . ^b^^bo^ 1*0, 

(c) Fast Flow (FF) ; 

= U/v^ ' = 1 Fig- <17c) 

P dp 

Space charge detaches from tip and 

^'fF^^'iF’ ^b^^bo increase again. 

Biblarz and Nelson (Ref. 16) stated that 
the maximum observed voltage ratio in the ordinary flow 
mode equaled but did not exceed the maximum reverse flow 
value of 1.13 for the 3.9 cm gap in their tests. This 
observation is consistent with the present model, since 
1.13 is the maximijm ratio for any flow mode and represents 
the entire physical gap. 

There still remains the strange behavior 
of the breakdown voltage curve for the short gap depicted 
in Figure (10) . Going back to Eq. (22) , recall that 
U/vdp = -^(1+c^), so it would be expected that at suffi- 
ciently short gaps and high curtents, convection effects 
would dominate the space charge at the tips, even at very 
' low velocities. This would lead to early (low velocity) 
separation of the space charge from the anode, resulting 
in an early rise in breakdown voltage probably to a near 
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since, as 



constant value at or below (^b'^^bo\max) 
decreased, d*— *-d. The data for the 1.6 cm gap appear to 
be a manifestation of this phenomenon. 

A numerical example may help to clarify 
the ideas outlined above. 

Let d=3 cm gap length 



V^^7800d+4000 = 27400 volts (18) 

Location of field point E=20,000 v/cm 
from Eq. (17) : 

^20K"”*^^ cm 

Effective gap length (due to space charge 
presence) : 



d' = d -5o = 3 - .45 = 2.55 cm 

Effective electric field (neglecting AV): 

E'=V. /d' = 27000/2.55 = 10,588 v/cm 
b o 

Assume Tf= 2(10) ^ sec, U= 76 m/sec = 

7600 cm/sec 

Anode extension due to convection: 

X = UT^ = 7600(2x10"^) = . 15 cm 
= d^-uT^ = 2. 55-. 15 = 2.40 cm 
and finally. 

Experimental values for V^/V^q the 76 m/sec velocity 
range fall in the . 94- . 96 range, with the current relatively 
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constant for all measurements. Above this velocity, con- 
vection at the tip, and/or the decreasing recombination 
distance (to be explained in next Section) takes over, and 
the effective gap length and begin to increase again. 

(5) Turbulence Effects on Breakdown Voltage . 
Figure (18) is a schematic of visually observed spatial 
characteristics of the tip glow patterns around the pins 
under no-flow, laminar flow, and turbulent flow conditions. 
It can be seen that the intense glow region of (a) is 
greatly extended by convection, (b) , and that the addition 
of turbulence tends to spread the glow over a much larger 
region. In fact, turbulence appears to shorten the maxi- 
mum extent of the glow under some conditions. It would 
seem probable that this spreading and reduction of the 
glow region would have an effect on the breakdown voltage. 
In fact, according to the model, it should increase 
due to physical shortening of 6 . This trend is well 
illustrated in Figure (10), where the curves for the 2,9 cm 
gap show a definite upward shift with increasing turbulence 
level. The sj^mbols, in order of increasing turbulence, 
are v , ° , , It should be noted that these symbols 

also indicate the order of increasing current, so again, 
increased current density could cause faster recombination 
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and smaller anode extension which would also cause an 
increase in V|^. Figures (19, 20) are photographs of the 
discharge section which illustrate the patterns of Figure 
(18). 

The next section analyzes in more detail 
the nature of flow and charged particle interactions in 
the cross -flow mode. 

(6) Cross-Flow Stability and Field Interactions . 
Figure (21) is a schematic representation of the space charge, 
electric and cross-flow field interactions as envisioned 
with the present model. The visual character of the dis- 
charge, as observed in many tests, appears to bear out the 
conclusions reached in the following analysis. 

Under no-flow conditions, the discharge 
gap experiences an apparent anode extension S^, due to 
space charge accumulation as in Figure (17a). As pointed 
out in the last Section, this apparent extension results in 
a reduction of effective gap length to d^=d- <S^. It was 
assiimed that the breakdown voltage was proportional to d^, 
so the maximum voltage ratio attainable in this configura- 
tion would be realized if the space charge could be removed 

completely, resulting in d*-*»d and V, /V, — ^d/(d- 6 ). 

o D bo o 
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Due to the nearly orthogonal orientation 
of the electric and flow fields in the cross-flow arrange- 
ment, these two effects operate nearly independently of 
each other to impart motion to any ionic space charge under 
their respective influences. For example, a positive ion 
formed at the tip of pin (a) in Figure (21) would experience 
electrostatic drift and convective motion as follows: 

X direction: 

Y direction: Vy=k"*’E(y) 

The important result is that the space charge location is 
strongly influenced by the magnitude of U. Taking the tip 
of pin (a) as the origin, the downstream displacement is: 

s = Ut (cm) 

U = gas velocity (cm/sec) 
t * = time (sec) 

For U = 30 m /sec (3000 cm/sec), and 
t = ^ = 2x10 ^ sec, the displacement x = 0.06 cm. For 
small tip radii, this amount of displacement can severely 
alter the local field structure, resulting in an as 5 mietric 
field distribution. Consider the region between the dis- 
placed space charge and the anode for example; this region 
would be characterized by a relatively weak field as the 
opposing positive point and space charge fields tend to 
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cancel each other. On the other hand, the region from 
which the space charge was removed will be at near the 
full xindisturbed point potential. The undistorted high 
field region on the upstream side of the pins would then 
be a source of continuous, stable corona, as electron 
avalanches feed into it and the space charge is quickly 
removed by convection. Visual observations of a cross - 
flow discharge give the impression that the stable, blue- 
violet glow, which characterizes region 1 of Figure (21), 
leans into the wind. This result at first seems contrary 
to intuition, but on closer inspection, it appears to be 
a natural result of the field interactions. 

On the other hand, while removal of space 
charge from region 1 enhances stability there, convection 
of that charge into region 2 promotes instability. 

It has been shown elsewhere (Ref. 7) that 
the presence of space charge enhances ionization in a uni- 
form gap if the reduced field (E/P) is less than one half 
the Stoletow constant B, which for air is 365 volts /cm- torr, 
and B/2=183 v/cm-torr. This is much higher than E/P in the 
non-uniform gap used in the present investigation, except 
in the immediate vicinity of the anode. It is evident then, 
that convection of space charge into region 2 will result 
in ionization instabilities there. 
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Manifestation of the instability described 
above takes the form of nxamerous arcs across the gap, com- 
mencing generally from the second and third rows of pins, 
and spreading to the first row of pins just prior to break- 
down. Even while heavy arcing is taking place from the 
first row of pins however, the arcs are confined to the 
region behind the blue-violet glow which marks the beginning 
of the discharge activity, thus demonstrating the destabi- 
lizing effect of the space charge. 

Since the power supply used in these tests 
is designed to "trip" on overvoltage and overcurrent (and 
it does faithfully in the ordinary flow mode) , the ntimerous 
arcs observed must be limited to fairly low currents, and 
it is not until applied voltage is raised to the value 
required to break down the relatively space-charge- free 
gap across the first row of pins that complete breakdown 
occurs, resulting in a large current pulse with resultant 
power supply shutdown. 

Actual breakdown streamer propagation 

7 8 

velocities have been reportedly measured in the 10^ to 10 
cm/sec range (Refs. 4, 10, 12), so that the breakdown 
streamer can travel the width of the gap (in the ranges 
tested) virtually unaffected by the gas velocity. For 
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example, in the longest cross-flow gap tested, d=3,2 cm, 

the time of transit of a primary breakdown streamer is 

t=d/v = 3,2/10^ = 3.2x10”^ sec. The displacement of an 

average molecule or ion by convection in this time inter- 

-3 

val, even at U= 150 m/sec is Ut = 5x10 cm. The result 
is that in the cross -flow mode, alteration of the field 
structure begins at low gas velocities by space charge 
convection, resulting in a stable upstream region of very 
low space charge, and a downstream region of instability 
caused by space charge distortion effects. Breakdown 
occurs across the nearly space-charge-free gap upstream 
at a value of = d/(d- 6^) > virtually unaffected by 

gas convection due to the high propagation velocity of the 
breakdown streamers . 

(7) Observed Effects of Charged Particle 
Density on Recombination Times . A very interesting phenom- 
enon was observed during the course of the cross -flow tests. 
While heavy arcing took place across all rows of pins, un- 
distorted by convection due to their high propagation 
velocity, a secondary pattern of Ixmiinous filamentary 
channels was also evident which appeared to be very much 
affected by the convective velocity (Fig. 21), Contrary 
to intuition again, however, the Ixmiinous secondary 
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streamers are convected a shorter distance at higher 
velocities. Their intensity also increased with increasing 
velocity. Figure (22) is a graph of recombination distance 
observations vs. velocity. 

In order to understand this behavior, it 
is necessary to focus on the exact mechanism which is 
responsible for these luminous streamers. In the streamer 
theory, the formation of a highly ionized channel from 
anode to cathode results in a highly conductive channel 
through which massive quantities of secondary electrons 
flow, having been released when the propagating streamer 
strikes the cathode. When the arc is extinguished, a 
column of space charge remains which is made up of positive 
and negative ions and electrons. When radiative recombina- 
tion takes place between particles of opposite charge, 
whatever the particular specific mechanisms, energy is 
released in the form of photons. The secondary streamers 
are believed to be caused by the massive recombination of 
the charged particles left in the aftermath of an arc. 

The standard approach to the recombination 
problem (Refs. 4, 7, 8', 13) is as follows: where charges 

of opposite sign reside together, they will be mutually 
attracted to each other, and so will come together, exchange 
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charge and neutralize each other, emitting radiation in 

the process. There are several varieties of recombination, 

all of which are discussed in detail in the references 

cited. It is sufficient for the present to say that the 

recombination rate will be at least proportional to the 

numbers of charged particles available, although the pro- 
0 

portionality varies with the process, so the change in the 

number of existing particles in a time (dt) will be given 

by an expression of the form: 

dn = - a n"^i dt (23) 

r 

where a ^ is a constant of proportionality. The nunnber of 
negative charges is not always equal to the positive charges, 
but n can be expressed as a fraction of n^, 
n” = cn"^ 

where c could also vary. Now, calling n"** = n, with the 

aid of the above, Eq, (23) can be written as 

dn = - a cn^dt 
r 



and 



dn/n^ = - a cdt (24) 

r 

which, when integrated yields the following: 

1/n - 1/n^ = a^ct 
showing that: 



t 



rec 



= 1/ a^c(l/n - 1/n^) 



(25) 
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Reference (13) states that recombination phenomena are 
usually recognized by the linear dependence of t on 1/n. 

But n is directly proportional to j , the current density: 
j = eCn'^v'*' + n"v”) = enCv"^ + v“c) (26) 

and also, it has been shown by other experimenters (Ref» 

14) and empirically from the present experiments that 

2 

j - U - n. Making use of this relation between n and U 
in Eq. (25) leads to the following: 
tree ~ 

Multiplying both sides by the gas velocity U gives the 
distance traveled during recombination in the arc-caused 
space charge column: 

x^ec “ ^^rec ~ 

which is the observed dependence. 

It is possible to do an order of magnitude 
analysis of by going back to Eq. (25). Assxmiing that 

c = 1, which is reasonable under these circumstances, and 
= 2(10) ^(cm^/sec) (Ref. 13). If it is assumed that 
the final ntamber of charged particles is = .In^ (90% recom- 
bination) , then 

tree ' a/2-10'^)(l/.ln^-l/n^) = (1/2-10'®) (9/n_,) (27) 

The relation between n^ and U can be 

o 

estimated by solving for the current in terms of n and U 
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as follows: 



j ^ en(k‘*‘E + k"E) (28) 

which reduces to 

j = en(k E) = env^ (29) 

since k^«k”. The drift velocity of electrons is difficult 
to estimate due to the variable field and severe space 
charge modification, but a reasonable approximation is the 
streamer propagation velocity - 10^ cm/sec (Ref, 12), 

This gives 

j = 1.6(10'^^) (10^)n^=1.6(10'^^)n^ (30) 

Now, from observed data, j(U) can be 

2 

calculated, assuming j - U . Taking a point from the data 
for atmospheric pressure and no turbulence; 
j = 15 ma = AU^ 

where typically U - 10^ cm/ sec. 

Then 

j = .015 amps = A(10^)^ 

A = 1.5(10”^^) (amp-sec^/cm^) 

and 

j = 1.5(10"^°)U^ (31) 

Now, equating (30) and (31), 

n^ = 90U^ (ions/cm^) (32) 
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This yields a value of = 10^-10^ ions/cm^, which is a 
reasonable value. 

Now, plugging this value of n^ back 

into Eq . ( 24) , 

= 1/2(10‘®)(9/90U^) and 

’'rec ° 5(10^) /U (cm) ( 33 ) 

4 

which yields for U = 10 cm/sec, 

o c 

Xrec = 5 cm 

It can be seen that the order of magnitude 
of the terms involved appears to approximately agree with 
the observed data. This tends to substantiate the basic 
assumptions regarding the secondary filaments as a recom- 
bination phenomenon. Figure (22) shows the comparison of 
✓ 

Eq. (33) with observed data. 

Assuming that the dependence illustrated 
in the preceding paragraphs holds for the "ordinary" flow 
mode, it can be seen that this idea can also be applied to 
the "virtual" anode concept in the following manner. Since 
Xrec a measure of the maximum extent of the space charge, 
it must be proportional to 6 , the anode extension, which 
is also a measure of space charge extent. This leads 
immediately to the following result: 
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Xrec - « ‘ a/U)C 

and 

''b * * > * ''bo<'*-=/U) (34) 

SO as U increases, increases. This effect has, of 
course, been observed as noted in the previous section at 
velocities above = 75m/sec. As a result of the above, it 
is difficult to ‘determine at this point which process, 
detachment of space charge from the anode or increasing 
recombination rate, or a combination of the two, is 
responsible for the upward swing of the ^^^^bo velocity 
curve in the ordinary flow mode (Fig. 23). 

(8) Reverse Flow Convection Effects . In this 
orientation (Fig. 9b), with no flow, the static anode ex- 
tension exists as in Figure (17a) . As velocity is increased, 
any positive ion located anywhere off the pin tip stagnation 
streamline will experience velocity components in the y 
direction due to the radial component of the anode electric 
field and in the negative and positive x directions due to 
convection and axial electric field components. The total 
velocity in vector form will be 

V = k'^E+U 

where the positive direction is from anode to cathode. 
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Due to the radially directed field lines 
near the pin tips, the resulting space charge pattern will 
tend toward the pattern of flow produced with a source 
flow immersed in a uniform flow field. The result is that 
velocity has a direct shortening effect on the anode ex- 
tension; but also note that the space charge tends to spread 
radially which reduces the charge density and therefore the 
field strength. The reduction of 5 and spreading tendency 
reduce the probability of streamer formation due to high 
charge concentrations, and thus reduce the current flow. 
Reduction in current means less charge, and so the glow 
collapses (this point will be covered in more detail in a 
later Section) , which quickly reduces the space charge 
distortion to the point that the breakdown voltage approaches 
the upper limit as 6 approaches zero (Fig. 24). It would 
be logical to assume that most, but probably not all, of 
the space charge anode extension would be destroyed or 

reduced below the critical level, and V, /V, ) „ = d/(d- s ). 

b bo max °o' 

It can be seen from Figure (14) for the 

reverse flow mode that this was the case. V, /V, ) for 

b bo max 

a 3.9 cm gap is predicted by the present model as 1.16 and 
the experimental value was 1.13 = ^b^^bo^max* 
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Werner (Ref. 15) reports on some work 
with a negative point-plane gap with velocity from the 
plane (anode) to the point (cathode) . The resulting space 
charge pattern and d 5 mamic processes should be similar to 
that described above, only with the negative ions playing 
the part of the positive ions in the former case. It 
might be expected then that the negative point would 
experience trends in the same direction as a positive 
point in regard to breakdown voltage and current changes 
with increasing velocity. It turns out that this is the 
case. Breakdown voltage increased with increasing velocity 
with a maximum = 1.08. Calculating for 

this case from Eq. (20) yields = 1.25, so the voltage ratio 
is less than the maximum. Reference (15) also reported no 
distortion in the glow pattern with flow while being ob- 
served under a microscope except for a reduction of size. 
This would be a logical result of the process described 
above for the positive point-to-plane since convection 
causes glow collapse, and an equilibrium glow region would 
be established which was much smaller than the no-flow 
glow. 

The apparent reduction of current at con- 
stant voltage in the reverse flow mode is of further 
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importance. Reference (16) reported no appreciable change 
in breakdown current in the reverse flow mode over the 
ordinary mode. This , imp lies a reduction in current at a 
given voltage, since increased appreciably; while 
Reference (13) definitely experienced lower currents at a 
given voltage with flow as the present analysis would 
expect. This fact will make it possible to determine the 
correct form for the functional dependence of the dis- 
charge current on velocity which will be undertaken in 
the next section. 

(9) Analytic Considerations in I = f(U) . While 

many investigations have derived expressions for current vs. 

2 

velocity of the form I = U (Ref. 14) , it can be shown that 
a more appropriate expression is I = (exp)^. 

In order to see this , it is necessary to 
examine some facts about the functional observed dependence. 
It is known that 

I “ f(U) Condition (1) 

and within practical limits of velocity, 

I - 0 Condition (2) 

Also, when U = 0, I = 1^ (the no-flow value). Condition (3) 
It has been shown also that when U <0 (reverse flow) the 
current tends to decrease. This implies that 
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dl/dU - 0 at U = 0 Condition (4) 

and from References (15, 16) , 

dl/dU =0 at U < 0 Condition (5) 

There are several f\mctions which could 
satisfy conditions (1 6e 2), but two of the simplest are: 



I 




I2 - (exp)*^ 



(35) 



To satisfy condition (3) , these can be written: 

11 . Au2 + 

1 2 - I (exp)™ 

^ o 



(36) 



Now, 



dlj^/dU = 2AU 
dl 2 /dU = I^C(exp) 



CU 



(a) 

(b) 



(37) 



It can be seen that Eq. (37a) does not satisfy condition 
(5), but (37b) does provided that C is a positive nxomber. 
Therefore, the equation which meets all the appropriate 



conditions can be written 
I = I^(exp)™ 



(38) 



In a later section it will be shown that, 
from a physical standpoint, the exponential form for 
velocity dependence evolves quite naturally, although in 
a slightly different form than Eq, (38). 
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B. GASDYNAMIC EFFECTS ON DISCHARGE CURRENT (NON-UNIFORM 

FIELD) 

As stated in the Introduction, a point- to -plane dis- 
charge can be divided into distinct regions due to the 
divergence of the applied electric field distribution. 

These are (1) an intense region of ionization near the 
point, (2) a space charge region fed from the ionization 
region by electrostatic repulsion of positive ions, and 
(3) the quasi-neutral region in which the applied electric 
field has dropped below the value necessary for appreciable 
ionization activity. This last region is disturbed more 
and more, however, as gap breakdown is approached by the 
intrusion of pre-breakdown streamers whose maximum lengths 
are proportional to applied voltage (Ref. 9). The actual 
length of any one streamer is a statistically varying 
quantity due to the nature of the dynamic processes in- 
volved. The amount of current associated with a specific 
streamer is proportional to its length, and current is 
also proportional to streamer frequency. Since both 
streamer length and frequency are proportional to applied 
voltage, the current must be nearly proportional to the 
sqxiare of the voltage. 

Several empirical expressions have been used to des- 
cribe the current -VO It age relationship in a point-to-plane 
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discharge under no-flow conditions. Reference (7), for 
example, gives the following empirical formula for a posi- 
tive point-to-plane: 

Io=CV(V-M) (39) 

where M is the starting potential of the corona current, 
and C is a fimction of the geometry and gas related 
variables. This formula gives the V dependence previously 
mentioned, and has been found to satisfactorily describe 
the I-V characteristic under no-flow conditions. 

Combining this expression with the current-velocity 
dependence developed in the preceding section, the 
equation is now 

I(U) * loCexp)*^” 

where is given by Eq. (39). can also be written as 
lQ=CkV(V-V^) which explicitly reveals the density dependence 
in the mobility k=kQ( p/ p ^) (Ref, 11) . This expression 
is also too limited, however, since it does not provide 
any means of introducing the effects of turbulence on dis- 
charge current. Since turbulence has proven to be a 
significant parameter, an equation to describe gasd 5 mamic 
effects on a discharge must be capable of reflecting 
turbulence as well as convection and density change effects. 
Such an expression can be found by going back to the basic 
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physical mechanism of current production in a gaseous 
discharge, the electron avalanche. 



In the special case of a uniform field with no space 
charge, the current (i) produced by an electron avalanche 
can be expressed as a function of the Townsend first 
ionization coefficient a in the following way (Refs. 4, 7, 



This equation indicates the exponential increase in current 
due to the nature of the electron avalanche process and is 
well described in the references cited. This expression 
becomes more complicated in a non-tmiform field where a 
being a function of the electric field, changes with x. 

It is then necessary to integrate a over the distance of 
avalanche travel: 



8 , 12 ): 



i = io(exp) 



(40) 




o 



(41) 



For purposes of the present discussion an effective 



can be defined over the avalanche travel distance: 




where x is the avalanche travel and 0^4 C -x, then 




( 42 ) 
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The subscript notation for a will be dropped for conven- 
ience, but it must be remembered that in a non-uniform 
field with space charge, a can vary in the space of a 
mean free path. 

A commonly used semi-empirical expression for a is 
the following (Ref, 4) : 

a/p = A(exp) ( 43 ) 
where A and B are constants of the specific gas xmder 
consideration. For air. Reference (4) gives 
A = 15 ionizations /cm- torr 
B = 365 volts/cm- torr 
with a reduced field range of validity of 
E/P = 100-800 volts /cm-torr. 

Of course, the current measured in a point -to -plane 
discharge circuit is the sxam of many small avalanches which 
form primarily near the point in the high field region; 
but streamers, which propagate far into the low field 
region are also built up from ntimerous secondary avalanches 
feeding into the tip of an advancing streamer. It seems 
apparent that any analysis of gasdynamic effects on dis- 
charge current, and hence power, must begin at this level 
to determine how the gasdynamics affect the electron 
avalanche . 
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1. Pure Convection Effects on Current 



Consider first the effects of pure convection. 

The model used in discussing breakdown voltage predicts 
that extension of the space charge region will initially 
occur as gas velocity is increased. This anode extension 
effectively moves the boundary of the high field region 
further from the anode, resulting in a longer average 
electron avalanche travel distance (Fig. 25). (Note that 
in discussing breakdown voltage, the apparent anode ex- 
tension was thought of as a reduction of effective gap 
length, and in discussing current effects, it is an in- 
crease in avalanche travel.) Recalling from the numerical 
example computed earlier the no -flow anode extension 6^ 
was 0.45 cm, and the apparent increase above this due to 
convection was 0.15 cm. An avalanche commencing at the 
beginning of the critical space charge region with con- 
vection would only travel 0.15 cm further than in the no 
convection case. This does not seem like much, but a 
simple calculation of avalanche current from Eq. (40) for 
both cases leads to surprising results. Using an average 

field strength of E =50,000 volts /cm, P=760 torr: 

a.v g 

a = AP(exp)^ ^^^^^^=44.4 ionizations /cm 



76 



1 

s 



V 



Ct 

and since = exp , for x = 0.45 cm (no flow); 

i/io>NF ' 4.77(10^) 
and for x = 0.45 + 0,15 (with flow): 
l/io)F = 3.71(10^^) 

and the increase of current for this single avalanche is 



i/io) 



= 80 



NF 



Although the numbers used above are only order of magnitude 
correct, the obvious conclusion is that even very small 
distortions of E-field structure by space charge convection 
can cause large changes in current. 

It is worthy of note at this point that the expon- 
ential current-velocity dependence derived in the previous 
section from purely empirical and functional considerations, 
turns out to be the same form the physical model would 
predict, that is an exponential increase due to increased 
avalanche travel in the high electric field region, 

2. Current-Voltage Characteristics 

As a further test of the space charge model pro- 
posed, it is instructive to look at other observed I-V 
characteristics. For example, the standard, no-flow I-V 
curve is known to be well described by the empirical 
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expression of the form of Eq, (39) . Now, in terms of 
physical processes, consider point-plane electrodes with 
increasing voltage. There will be no ionization activity 
mtil the point reaches the potential required for ioniza- 
tion of the surrounding gas molecules, which gives rise 
to the "inception” voltage (M or V ) , at which the corona 
current starts as small weak pulses caused by the initial 
formation of the electron avalanches. As voltage continues 
to rise, the high field front of the electric field moves 
out into the gap, causing the average distance at which 
avalanches can form to move with it. As already demon- 
strated, this results in a rapid rise in current, which, 

from a physical standpoint, is better described by an 

2 

exponential function rather than a U type of empirical 
fit. 

A simplified representation can be obtained for 
the total current by observing the general form of the 
result of a summation of several individual avalanche 
currents. This will be useful in later developments and 
for a general trend analysis to further test the validity 
of the proposed model. 

Asstime that the total number of avalanches per 
unit time is N. Then the total current is 
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1 



= + i2 + . . . 



I 



N 



= E 

j'l 






Since most activity takes place in the same small region 
of the pin tip, let the following average quantities be 
defined: 



N 

X = (1/N) E X, 
j=l ^ 



N 

a = (1/N) a . 

j=l J 



= (1/N) 



N 



E 

j=i 




Sxmming up all the avalanche contributions, the 
total current is 

I^=i^ + i 2 + ... = iQ(exp) + i^(exp) + ... 



_ _ , .ax _ 

Io=Ci(exp) - 

where C, = i N and C« is a function of the electrode 
1 o Z 

geometry, since no exponential increase in current can 
take place until the electric field strength at the electrode 
surface at least reaches the ionization strength of the gas 
surrounding it. Discounting normal saturation currents 
due to free electrons which are always present in the air. 
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C 2 can be evaluated by the fact that = 0 at the inception 
of the corona current: 

Iq - 0 = C^(exp) “ - C 2 



SO 



C 2 = C^(exp) 



at corona ignition 



but at the point when ionizing activities first commence, 

X = 0, and 

C2 “ C^(exp) ~ ^1 ~ ^ 
and, therefore, 

= C(exp“ ^ -1) (44) 



Equation (44) is the basic physical result for non- 
imiform fields under no-flow conditions as derived from the 
space charge model. The associated constant is again a 
function of gas type, density and electrode geometry. It 
is interesting to note that this is nothing more than the 
expression for current between parallel plates in a uniform 
field, written to accommodate variable geometry. 

Closer inspection of the constant C = i^N shows 
that it is still a function of streamer length L = f(V) 
and frequency N = f = V, so can be rewritten as CV"^ and 
then Eq. (44) can be written as 

= Cv2(exp -1) 
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or, taking one V inside, 

a X 

= CV(Vexp - V) 

and from the condition that = 0 and x = 0 at V = V^, 
this results in 

a X 

= CV(Vexp -V^) 

This is nearly identical with the previous empirical forms 
for the no-flow condition, except that now the exponential 
expression introduced permits easy interpretation of gas- 
dynamic effects on the discharge current. 

The X in Eq. (44) is proportional to 6 , so it will 
be less confusing to rewrite the current equation in the 
following form: 

= CV(Vexp^ (45) 

remembering that the subscript 'o' represents no-flow 
conditions. Equation (44) is plotted in Figure (26), 

It is now easy to see the effects of increasing 
velocity on the current, since in the ordinary mode, 
increases with flow. This would increase the exponent 
which has the effect of shifting the I-V characteristics 
up and to the left as in Figures (27 a&b) . These curves 
are typical of the characteristic shifts which have been 
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observed in the present experimental work. The general 
form of 6 in terms of the other variables is 
6 = 6q + UT - f(u) = c— + UT -f (u) , where f(u) is the 

result of space charge "blowout" and recombination effects. 

Looking again at the reverse flow case discussed 
earlier with the aid of Eq. (45), it is now easy to see 
how a low current equilibrium situation is established. 

As soon as there is even a slight distortion (shortening) 
of the space charge region, the current decreases along an 
exponential curve, which results in less charge production, 
further decrease in current, etc., until equilibrium is 
again established, when the space charge region shrinks 
far enough into the high field region of the anode that 
it is no longer affected by convection. Figure (28) 
summarizes the effects of velocity on 6 and the breakdown 
voltage characteristic curves, 

3, Diffuser/Pensity Effects on Current 

The effects of reduced density on the current are 
also easily deduced from Eq. (45), since as shown earlier, 

5 increases with a density reduction, but also a (p) 
increases with decreasing density over the ranges of 
reduced pressure discussed so far. (There is an optimum 
ionization efficiency, n a/E, which occurs for air at an 
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E/p = 365 volts/cm-torr = B, which will be discussed below 
(Ref, 12)) . The result of this dual exponential effect is 
a rapid rise in current with decreasing density, until the 
reduced pressure in the ionization region drops below the 
optimum. Beyond that point, the effects of density reduc- 
tion on the current could actually reverse if the ionization 
efficiency drops too far (Fig. 29) . 

To arrive at an estimate for optimum pressure for 
best ionization efficiency, it is possible to write 
E/P)^p^= 365 v/cm-torr = B 

and using E = 200,000 volts /cm (value at pin for a 3 cm 
gap, h}TJerboloidal point-to-plane) , 

Popt= E p/B = 200,000/365 = 547 torr 

The value computed from the data of Reference (18) 
in Appendix (b) for a 3 cm gap is = 480 torr, which is not 
too far off, considering the approximation inherent in the 
point model used. (For example, if the point radius were 
assumed to equal the wire radius of = 0.08 cm, then 
Popt”375 torr) , Since the discharge in the present 
investigation was operated at a pressure well above the 
optimum value for most efficient ionization, it would be 
expected that any density reduction produced in the test 
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section by the addition of a diffuser would enhance ioniza- 
tion significantly, which it does (Fig. 30). 

4 . Turbulence Effects 

While the study of turbulence is, in itself, a 
complicated field, the basic interactive effects of turbu- 
lence on an electric discharge can be qualitatively des- 
cribed with only the most basic characteristics of a 
turbulent flow field; those of fluctuating, time -dependent 
variations in gas density and velocity* 

The turbulence generated in the discharge section 

o 

in the present work ranged in frequency from = 100-20,000 Hz 
(Ref. 3). This frequency range yields characteristic times 
= (1/f) of the order of 2(10"^) -5(10”^) sec. These 
times are long compared with the characteristic times 
associated with charged particle kinetic processes (Table I) . 
The result is that the charged particle properties adjust 
in a quasi- steady fashion to fluctuations in gas properties. 
Thus, charged particle kinetics and energy transfer pro- 
cesses couple to fluctuations in the neutral gas primarily 
through changes in gas temperature and density (Ref. 19). 

For example, consider a point in a turbulent dis- 
charge region; gas density at this location will vary 
periodically around a mean value. Although the time 
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average of the fluctuations at the point is zero, and the 
point experiences as much time at elevated pressures as 
it does at reduced pressures, the net effect of these 
fluctuations on discharge current are very significant. 
This is due to the exponential dependence of the local 
value of the ionization coefficient on the local reduced 
pressure (E/P). As long as the discharge (E/P)> B, 
pressure fluctuations will significantly increase the aver 
age effective ionization coefficient. 

EXAMPLE: P = 760 torr = mean pressure 

p = 50 torr = fluctuating component 
ct = AP(exp) 

and using a simple expression for current for illustration 
I = 

Computing ionization coefficients for E = 50,000 v/cm: 



“760 


= 44.4 


Mean 


“810 


= 32.9 


High P fluctuation 


“710 


= 59.8 


Low P fluctuation 



Current computations for = 1, x = 0.5 cm: 

I = (exp)^^*'^^’^)= 4.38(10^) Mean 

I = (exp)^^*^^*^^= 1.39(10^) High P 

I = (exp)^^'^^*^^= 9.67(10^2) Low P 

lavs' Uhl + ' ^.83(10^2) 
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The amplification of the current over the mean due to a 
fluctuation in pressure is then: 

I _=4. 83(10^^) /4. 38(10^) = 1103 

avg mean ' ' ' 

Obviously, the effects of turbulent fluctuations in pressure 
or density can cause substantial increases in current due 
to the quasi-steady coupling of electron kinetics to temporal 
variations of the neutral gas properties. 

Local increases in electron/current density, however, 
can have a great destabilizing effect on the discharge, 
since increased local current density causes increased 
joule heating (j*E) of the gas. This heating results in 
local density reduction, increased conductivity and further 
increases in current density which may result in a runaway 
condition leading to complete breakdown of the gap (Ref. 20) . 

Fortunately, it is the basic nature of turbulence 
that the strongest density fluctuations are accompanied by 
the highest fluctuating velocity components (Ref. 19) . 

These high velocity components tend to disperse concentra- 
tions of space charge, dissipate local hot spots and pro- 
mote stability. The role of turbulence, then, is two-fold: 
on the one hand, it enhances ionization, while on the other 
hand, it jointly acts to prevent an unstable, rtinaway 
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condition from developing due to formation of local hot 
spots in the gas. 



Before leaving the subject of turbulence, it is 



worthwhile to do a general frequency analysis to determine 
how a given turbulent frequency spectrum will affect a 
discharge. It was pointed out in Reference (16) that 
turbulence spectra which were heavily weighted in the low 
frequency regions gave the best improvement in power enhance- 
ment to the discharge. 



effects of turbulence more explicitly on the ionization 
coefficient a . In a turbulent field, a = a(t) , and if 
f ( <ij t) is a periodic fluctuation of gas density, then the 
time average of a(t) over one cycle will be: 



but since o) is an independent parameter of the problem, 
it may be included as follows: 



This shows that a is inversely proportional to the turbu- 
lent frequency components. It should be possible to 
develop a weighting factor for a given turbulence spectrum 



That this should be so can be seen by analyzing the 





( o> t) d( a> t) = — 



(46) 
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m 



by siinming the contributions of the various frequency 
components to the ionization coefficient. It must be 
realized, however, if an attempt is made to quantify the 
effects of a given turbulence spectrum, that inhomogenie ties 
in the turbulence pattern can lead to premature breakdown 
through a relatively inactive "dead" region of the gap 
where instabilities may develop uninhibited. 

The following is a summary of the gasdynamic effects 
on discharge current observed in the present experimental 
work and interpreted in terms of the present model: 

a. No -Flow 

I = CV (Vexp*^l“'^°-V-) 

O ' t' \ c 

Current commences at an inception voltage 
and rises on a modified exponential curve until breakdown 
occurs (Fig. 26) . 

b. Flow (ordinary mode) 

I(U) = CV(Vexp^“'^ ) 

c 

The anode extension 6 increases due to con- 
vection of space charge, shifting the characteristic curve 
up and left, while a also increases due to high space 
charge in the gap (Fig. 27b) , 
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c. Reduced Density (diffuser on) 

6 and a increase due to reduced density, and 
6 also increases due to the convection effect. The result 
is a large shift of characteristics up and left. The re- 
duced density also reduces the breakdown voltage signifi- 
cantly, and may result in a reduction of maximum power 
(I^V^) (Fig. 30). 

d. Turbulence (diffuser off) 

Substantial improvements in a can be realized 
with the proper t}rpe of turbulence, and stability is con- 
currently enhanced. Mean velocity also causes 6 exten- 
sion due to convection, and a large increase in a 6 
results, with the accompanying shift up and left of the 
I-V curve (Fig. 27a). 

e. Turbulence and Density Reduction (diffuser on) 

The combined effects of (c) and (d) above 

result in maximum I-V shift and high power capacity due to 
stability enhancement by turbulence. 

Figures (19, 20) are a sequence of pictures of 

1 

the discharge section and the discharge operating under 
some of the conditions described above, and Figure (18) is 
a schematic of the pin tip regions as they appeared during 
visual observation. One additional feature of turbulence 
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which is evident from these figures is that in addition to 
homogenizing the discharge region the turbulence creates a 
much larger ionization region at the anode, which is the 
major source of the additional current. 

5 . Maximum Current/Discharge Stability 

Figure (31) shows how the breakdown current varied 
in the present investigations as a function of mass flow 
rate through the discharge region for the various con- 
figurations already mentioned. 

The curve for pure convection (no diffuser^ no 
turbulence grid) shows increasing stability with increasing 
mass flow rate. This stability increase is probably due, 
at least in part, to convection of heat out of the high 
current density regions near the pin tips where growing 
thermal instabilities may develop under no-flow conditions. 
The data indicate that a limit is apparently being approached 
in this mode. This would be expected in a relatively smooth 
laminar type of flow, since the joule heat developed at a 
pin tip would flow toward the cathode with relatively slow 
radial diffusion, resulting in a coltmin of heated gas be- 
tween anode and cathode, increasing the conductivity and 
breakdown probability of the gap. 
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The addition of the diffuser raises the curve, 
thus indicating a stabilizing effect associated with re- 
duced density. This stabilization can most probably be 
traced to the increased diffusion rates of charged parti- 
cles, reducing high density charge build-up, and increased 
thermal conductivity, which also tends to reduce local 
heating. The problem of laminar flow is still present, 
however, and it can be seen that, once again, a current 
or power limit is being approached as mass flow rate 
increases . 

The addition of the turbulence generating screen 
in either of the above configurations demonstrates a very 
marked increase in stability, with no apparent saturation 
limit being approached within the limits of the power 
supply capabilities used in the present investigations. 

It thus appears that the fluctuating, large scale, high 
intensity velocity components introduced are very effective 
in enhancing discharge stability. This enhancement is 
presumed to act through intense mixing of ionic space 
charges with surrounding areas, and rapid dissipation of 
heat from local regions of the discharge gap. The result 
is a smoothing out of the electric field distortion caused 
by space charge, and a rapid redistribution of heat 
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throughout the entire discharge volume, whereas in the 
laminar flow cases, the heat was confined to relatively 
small coliimnar voltimes between anode and cathode. 

The most effective screen design turned out to be 
a redesigned "number 9" screen of Reference (6) « Reference 
(3) gives detailed engineering data on all equipment used 
in the present work. One of the significant features of 
the redesign was the incorporation of large boxmdary layer 
bleed areas around the discharge volume (Fig. 2). Flow 
velocities measured through this boundary layer region 
were 2-3 times the average discharge velocity with the 
turbulent screen in place. It is felt that this high 
velocity region serves as a heat sink for the discharge 
and greatly enhances the thermal stability by very rapid 
removal of waste heat generated in the discharge region. 

Reference (3) should also be consulted for an 
interesting approach to power input optimization for the 
type of flow device used in this work. 
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V. CONCLUSIONS 



It appears that the space charge convection model 
presented here is a valid approach to electric discharge 
and gasdynamic interaction analysis in highly non-uniform 
fields. The model has also been shown to hold as an 
approximation to discharge breakdown properties in a 
iiniform field, and so may not be limited to non-uniform 
geometries where it appears to describe many observable 
effects concerning gasdjmamic interactions. 

Gasdynamic interactions promote stability of a dis- 
charge in several ways, through pure convection, density 
reduction, and turbulence. It appears that the basic 
mechanism of all three of these is through thermal stabi- 
lization by heat removal from potentially developing ’’hot 
spots” in the gas. 

High intensity, low frequency turbulence is best for 

« 

stabilization. This fact may permit quantitative analysis 
of individual turbulence spectra for use in optimization 
studies. 

The addition of a diffuser and turbulence generating 
screens results in much higher power input at lower mass 
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flow rates, and therefore is an important economic 
consideration in the construction of larger devices of 
this type. 
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VI . RECOMMENDATIONS 



A. The major recommendation resulting from this study 
is that, initially, more comprehensive testing should be 
done on all aspects of the gasd 3 mamic/flow interactions 
reported in order to establish accurate empirical inputs 
required for optimization studies. Concurrently, attempts 
should be made to derive accurate relations among the 
parameters discussed using the Poisson Equation and 
appropriate flux equations. 

B. The stability problem is probably the most critical 
and complex, but since it is the key limiting performance 
parameter for high pressure discharges, an accurate means 
of quantification is necessary which properly reflects the 
influence of gas d)niamics. It is hoped that the relatively 
simple model presented here will aid in visualization of 
the stability problem and the gasd}mamic factors which 
affect it. In regard to discharge stability, it is 
specifically recommended that the following tests be 
performed: 

1. Current distribution studies should be performed 
for correlation with temperature profiles and breakdown 
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location statistics, and turbulence and velocity profile 
measurements . 

2. Oscillographic studies of current and voltage wave 
form characteristics should be attempted in the various 
configurations tested, and at various gap lengths to deter 
mine the exact type of breakdown which occurs . 

3. Small diameter wire should be tried as the anode, 
with the wires mounted parallel to the cathode. Wire 
electrodes resulted in a much more uniform discharge in 
recent MHD laser cavity experiments (Ref. 22). 

C. Since the present work was limited by power supply 
capabilities, a larger power supply with much greater 
power output is necessary to extend the investigation of 
optimum gap length, and maximum power input limits, 

D. An attempt should be made to determine more specif- 
ically the mechanism responsible for the increase in at 
high velocities in the ordinary flow mode, since small 
increases in at that point in the curve result in a 
large increase in due to the slope of the I-V 
characteristic. 

E. Extension of this work to a lasing mixture should 
be attempted when a better power supply is available. It 
is believed that a lasing mixture may prove more thermally 
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I 



stable, since most of the input energy goes into the 
vibrational molecular states instead of translational 



gas heating, 
coupled into 



and therefore even more power can be stably 
the gas. 
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1,50-xiets (kv) 




Figure 1. Sample Electric Field Distribution Curves for 
Hyperboloidal Point-Plane and Concentric Spheres. 
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(d) Boundary Layer Slot Screen.. Used with test section of 
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Figure 4. Power Supplies 
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Figure 5. Variable Speed Transmission 
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Figure 6. Electric Field and Voltage Distribution between 
Parallel Plates. Subscripts: (1) No Space Charge 

(2) Moderate Space Charge (3) High Space Charge 
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GAP LENGTH d(cm) 

Figure 7. Measured and Calculated Values of E|^ for Parallel Plates 
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Figure 8. Predicted Model Breakdown Voltage Ratio 
Density Dependence vs. Actual Data. 
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(a) Ordinary Mode 
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(b) Reverse Flow 
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(c) Cross-Flow 



Figure 9. Various Orientations of Electrodes and 
Turbulence Grids Tested. 
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Test Section 1 




Figure 10. Voltage Ratio and Current vs. Velocity 
in the Ordinary Flow Mode. 
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Figure 11. Voltage Ratio vs. Velocity in Referse 
and Cross -Flow Modes. 
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Figure 12. Electric Field Distribution in a Point-Plane 
Gap. Curves: (1) Point Field Only (2) Space Charge 

Only (3) Total Field 
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Figure 13. Breakdown and Corona Initiation Voltage vs. Gap Length 
for No-Flow at Atmospheric Pressure. 
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Figure 14. Summary of Maximum Breakdown Voltage Ratio vs. Gap Length 
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Figure 15. Calculated Values of the Slip Parameter at 
the Field Point E = 20,000 volts/cm. 



112 



i'i 





Figure 16. Calculated Values of Slip Parameter at Anode 
Surface. Note that Sp = f(d) . 
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Figure 17. Various Stages of Anode Extension and "Blowout” 
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(0) NO FLOW 




(b)FLOW, NO TURBULENCE 




Figure 18. Visual Character of Various Glow Modes 
without Diffuser 
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Figure 19. Photograph of Discharge. No Flow 

(10 sec exposure). Pins on Left. Bright region 
on bottom is an arc which took place during 
timed exposure. Note small spherical tip glow 
region. 
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Figure 20. Pins on Left. Ordinary flow from left to 
right (1 sec exposures) . Note elongation of glow 
region in (a) due to convection (u - 150 m/sec), and 
very homogeneous spreading of glow in (b) due to high 
turbulence. 
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Figure 21. Schematic Representation of Cross-Flow Current Distribution (Observations). 
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Calculated and Measured Recombination 
in the Cross -Flow Mode. 



150 



119 






I 

h 




o 



> 

o 

o 

UJ 

> 






o 



o 

M-l 



0) 

T) 

O 

S 

>> 

rO 

T3 

0) 

U 

o 

•H 

(U 

S-J 

pL. 

CO 

C\J 

CO 

a 

•H 



CO 

•H 

5-1 

0) 

4J 

O 

ccJ 

5 -» 

CCJ 

4i: 

o 

o 

PO 

CCJ 

4J 



0 

> 

a 

1 

X) 

Cd 

o 

5-^ 

PQ 



cd 

5 -» 

(D 

C 

(U 

O 



cn 

(N 

CD 

5-1 

D 

W) 

•H 



120 



Ordinary Flow Mode. 
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X 

t ® ^ Electric Field 

(1) No Space Charge 

(2) With Space Charge, No Flow 

(3) With Space Charge, and Flow 
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(b) Schematic of Flow Effect on Space Charge 



Figure 25. Convection Effects on Average Electron 
Avalanche Travel Distance 
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Figure 26. I-V Characteristic Velocity Dependence from Eq. (44) 
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Figure 27. Convection and Velocity Effects on I-V 
Characteristics. (a) 3 cm gap, with Turbulence, 
(b) 2.9 cm gap, no Turbulence. 
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Figure 28. Summary of Flow Effects on Breakdown Voltage with Representation of 
Pertinent Model Parameters . 
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Figure 29. 
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(Modified Scr #9 Used for Turb) 
O Diffuser On, With Turb 

□ Diffuser Off, With Turb 
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Figure 31. Breakdown Current vs. Discharge Mass Flow Rate for Various Configurations. 



APPENDIX A 



VARIABLE VOLTAGE RATE TESTS 

As stated in the Introduction, this series of tests was 
vmdertaken in an attempt to determine whether any type of 
resonance condition might exist between the rates of voltage 
increase used and any of the d}niamics of the discharge or 
gas parameters. The fact that the characteristic times for 
current pulsing, rates of voltage increase and turbulence 
all overlap in the ranges tested was the primary motivating 
factor. 

The breakdown voltage was found to decrease in a linear 
fashion with increasing voltage rate ( y ) in the ordinary 
and cross-flow modes (Fig. A-1) . At higher velocities in 
the cross-flow mode, the breakdown voltage increased at the 
very low values of y as would be expected from the space- 
charge -convection model. The slope of the curves remains 
basically unchanged, however. 

In the ordinary flow mode, the effect of the flow 

definitely increased the slope of the V, curve, but V still 

b b 

decreased linearly with increasing voltage rate. 

Under no-flow conditions, using the Townsend type of 
formula for the relation between current and voltage, it 
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can be shown that the slope of these vs. 7 character- 
istics should be negative. Assxjme is a linear fxmction 
of y , and 

where the 'o' indicates the value at >^ = 0. The Townsend 
current equation for breakdown is 

Ib ' - V,,) ■ 

and breakdown power is P|^ = 

Pb ' ' '^(''bo+A )^-C(Vb„+A 

and 

dPj,/d = 3C(V^;,+A )2 a-2C(V^^+A )AV^ 

Now, assuming has a maximirai or minimum, letting the 
above = 0: 

max/min ' <2V,/3 - /A 

Since is always positive and 2V^/3, then A, the 

slope of the curve, must be negative as observed. 

The reason for this decrease in V, and I. is not known 

D b 

and was not pursued due to the uninteresting trends observed. 
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Cross and Ordinary Flow Modes. 



APPENDIX B 



NO -FLOW EMPIRICAL GAP LENGTH/DENSITY 
OPTIMIZATION AND COMPARISON WITH OBSERVED DATA 

Prior to coinmencement of the testing phase of this work, 
data was gathered from a nvimber of sources on the break- 
down properties of a point-to-plane gap. Empirical equa- 
tions for ) and ) were developed and used in 

a computer optimization study for this electrode configuration. 

The results of this study indicated that, under no-flow 
conditions, maximimi power /cm input could be obtained at 
d = 4 cm and a density ratio of 0.65 = ( / ^) . This 

corresponds to a pressure of = 420 torr (Fig. B-1) . Figures 
B-2 & B-3 show the computer generated and I^ character- 
istics plotted against density ratio. These figures indi- 
cate that the 4 cm gap was best at atmospheric pressure 
too, while data from the present experiments yielded an 
optimvim gap of = 1.6 cm under no -flow conditions at atmos- 
pheric pressure (Fig. B-4) for a multiple pin electrode. 

The equations used in the optimization study are listed 
below for reference: 

Breakdown Voltage: 

V,( ,d) = 51( / / )(l-eKp’^)+B 

D O o 
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where 



T = d^*^/3 

B = 20(d/d^-l)( / / ) 

d = 4 cm 
o 

Breakdown Current: 

Ib( ,d,Vb> = .09C(( / - A) 

A = 590(d/dQ-l)( / -.04)'^^(1.2- / ) 

E = 1.05 + ( / - .05)*^^ 

C = 1 for d 4 cm 

C = .lld^-^ for d 4 cm 

A series of tests in the ordinary flow mode was con- 
ducted using the turbulence screen number '’9" (Fig. 2) of 
Reference (6) at various velocities and gap lengths in an 
effort to determine the optimum gap length at atmospheric 
pressure in the present apparatus. Figures B-5 through B-8 
show the results of these tests. 

In Figure B-5 the typical vs. velocity trend is evi- 
dent, showing an initial decrease of and a later increase, 
except for the short gap-high current case of the 1.6 cm 
gap which was explained earlier. 

Figure B-6 is interesting in that it shows reaching 
a maxim;jm value of = 17 ma for the three best gaps. This 
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may be due to the fact that, with no boundary slots as in 
later tests, the blockage of this screen is = 60%, and maxi- 
mtim velocity is limited to = 120 m/sec through the discharge 
section. This may result in the section becoming thermally 
unstable due to inefficient heat removal. 

The maximtim specific power is shown in Figure B-7. Here 
it can be seen that the short gap permits highest specific 
power. This would be expected, since the shorter the gap, 
the lower the maximum temperature rise in the gas at a 
given velocity (Ref. 5) . 

Figure B-8 shows the breakdown current variation with 
gap length for various constant velocity profiles from the 
data of the previous figures. It is interesting that at 
high velocity and high turbulence, the 4 cm gap current 
increases faster than at the other gap lengths. 

It must also be noted that it was later discovered that 
the power supply used in these tests was limited to 18-19 ma 
by internal problems, so part of the apparent 17 ma limit 
in the discharge illustrated in Figure B-6 may be attributable 
to a power supply limitation rather than a discharge stability 
limit. 
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DENSITY RATIO 

Figure B-1. Breakdown Power/cm vs. Density Ratio for Single Point-to-Plane 
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Figure B-2. Breakdown Voltage vs. Density Ratio for Single Point-Plane 
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Figure B-3. Breakdown Current vs. Density Ratio for Single Point-Plane. 
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Figure B-4. Breakdown Power/cm vs. Gap Length (No Flow) for Multi-pin to Plane 
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U(m/s) 

Figure B-5. Breakdown Voltage vs. Velocity for Multi-Pin Configuration with 
Turbulence Screen #9. (Various Gaps) 
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Figure B-6. Breakdown Current vs. Velocity for Multi-Pin Configuration with 
Turbulence Screen #9 
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with Turbulence Screen #9 
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GAP LENGTH (cm) 

Figure B-8. Breakdown Current vs. Gap Length for Multi-Pin Configuration with 
Turbulence Screen #9 
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